Energy Transduction of Isothermal Ratchets: Generic Aspects and Specific Examples 

Close to and Far from Equilibrium 
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We study the energetics of isothermal ratchets which are driven by a chemical reaction between 
two states and operate in contact with a single heat bath of constant temperature. We discuss 
generic aspects of energy transduction such as Onsager relations in the linear response regime as 
well as the efficiency and dissipation close to and far from equilibrium. In the linear response regime 
where the system operates reversibly the efficiency is in general nonzero. Studying the properties 
for specific examples of energy landscapes and transitions, we observe in the linear response regime 
that the efficiency can have a maximum as a function of temperature. Far from equilibrium in the 
fully irreversible regime, we find a maximum of the efficiency with values larger than in the linear 
regime for an optimal choice of the chemical driving force. We show that corresponding efficiencies 
can be of the order of 50%. A simple analytic argument allows us to estimate the efficiency in this 
irreversible regime for small external forces. 

PACS Numbers: 87.f0.+e, 05.40.-a 
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I. INTRODUCTION 

Biological systems provide an important motivation to 
study the physics of active processes which on a molecular 
scale are able to transduce chemical energy into mechan- 
ical work and motion. Important examples are motor 
proteins and enzymes which move actively along DNA 
0. The properties of such systems differ in several re- 
spects from macroscopic machines and heat engines: (i) 
active phenomena occur on a molecular scale in a very 
viscous environment with overdamped dynamics, motion 
is thus stochastic and obeys only on average the first and 
second laws of thermodynamics; (ii) these systems are 
isothermal and operate strictly at constant temperature 
as they are in intimate contact with a thermal bath. In 
recent years, a number of theoretical approaches to de- 
scribe this class of systems have been developed [^-^ . 

In order to discuss the energy transduction of such sys- 
tems, the concepts which have been developed for macro- 
scopic motors have to be applied with some care. Re- 
cently, there has been a growing interest in the energetics 
of Brownian motors p-pl| . It is the aim of this article to 
discuss generic aspects of energy transduction of Brown- 
ian motors driven by a chemical reaction and to provide 
several specific examples which reveal new and interest- 
ing properties. 

The two-state models which we use § represent a use- 
ful paradigm for the description of energy transduction 
of isothermal motors in the overdamped regime. They 
are motivated by cytoskeletal motor proteins which move 
along polar and periodic filaments. Coupling a two state 
model to a chemical reaction, which induces transitions 
between the two states of the motor, leads to motion and 
force-generation if the chemical potential difference A/i 
between the fuel and its reaction products is nonzero and 



if the system has a polar symmetry. Assuming that the 
chemical reservoirs coupled to a single motor are macro- 
scopic in size, this chemical potential difference can be 
defined even under out-of equilibrium conditions since in 
this limit the reaction driving the microscopic motor af- 
fects the reservoir only weakly. Using A/x as the relevant 
control parameter, the consumed chemical free energy by 
the active process is well defined. This leads to a simple 
definition of efficiency ry as the ratio of the mechanical 
work performed and the consumed chemical free energy. 
We find three important results: 

• The efficiency calculated for these models can be 
maximized far from equilibrium. 

• Close to thermal equilibrium there exists a linear 
response regime which is important because of its 
universal features. We demonstrate that the de- 
pendence of the efficiency in this regime on tem- 
perature is strongly model dependent and can be 
non-monotonous in which case thermal fluctuations 
are essential for an efficient energy transduction. 

• The efficiencies vanish at stalling conditions ( zero 
average velocity) except in a singular limit where 
they reach the ideal value 77 = 1 . 

The outline of our paper is as follows. In section II, 
we discuss generic aspects which are completely inde- 
pendent of the model chosen. We define the efficiency 
and identify the generalized currents and forces which 
allow us to write a linear response theory. We discuss 
the generic features of efficiency in this regime, in partic- 
ular the maximal efficiency under reversible conditions 
and the efficiency at stalling conditions. In Section III 
we choose an explicit realization of the transport equa- 
tions where the motor is described as a two-state model 
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which is coupled to a chemical reaction and we identify 
the energy fluxes in the system. Section IV discusses the 
energy transduction properties for specifically chosen ex- 
amples. We show that efficiency is typically optimized in 
the irreversible regime and give examples for the temper- 
ature dependence of 77 when the system operates in the 
linear response regime. In our concluding remarks, we 
relate our results to biological motors and discuss alter- 
native definitions of efficiency which have been used in 
the literature. 



Because of energy conservation, the amount of energy 
dissipated per unit time therefore reads: 

n = ,Utv + rAfi . (7) 

From the second law of thermodynamics it follows that 
n must always be positive. 

B. Linear response theory 



II. ISOTHERMAL RATCHETS: GENERIC 
ASPECTS 

A. Force, velocity and efficiency 

Motivated by linear biological motor proteins which 
move along a linear filament, we will consider chemically 
driven systems which can induce motion along a one- 
dimensional track. The energy source is the difference of 
the chemical potentials A/i of fuel and products. Being 
motivated by biological motors, we use the hydrolysis 
ATP ^ ADP -I- P as example |1. We define 



(1) 



where fiA and /ip are the chemical potentials of ATP 
and ADP-I-P, respectively. In order to perform useful 
mechanical work, the system has to move against an ex- 
ternal force /ext applied parallel to the track. In addition 
to the two generalized forces A/i and /cxt acting on the 
system, we can define two generalized velocities: (i) the 
average velocity of motion v of the motor along the track; 
and (ii) the chemical reaction rate r defining the average 
number of ATP molecules consumed per unit time. The 
motor can thus be characterized by the equations of state 



r = ^(/oxt, A//) 



(2) 
(3) 



which describe the velocities of the system as a function 
of the generalized forces . The mechanical work per- 
formed per unit time against the external force is given 
by 



W = fcxtV 



(4) 



The amount of chemical energy consumed per unit time 
is 



Q ^ rA/i 



(5) 



For a system which performs mechanical work, i.e. 
foxtv < 0, we can define the (mechanical) energy trans- 
duction efficiency as psl 



V = 



rA/i 



(6) 



Close to thermal equilibrium, i.e. for small forces 
/cxt ^ T/l and A/i ^ T, where I is a typical length 
scale of the motor and T is the temperature measured in 
units of ks, we can expand Eq. (H) to linear order: 



All /cxt 
A2l/oxt 



Ai2A/i 

A22A/i 



(8) 



The matrix of linear response coefficients has the fol- 
lowing physical meaning: An is a mobility giving the 
response of the velocity to the applied force. A22 plays a 
similar role for fuel consumption. It describes the 'chem- 
ical admittance' or the response of the chemical reaction 
rate r to the chemical force A/i. The coefficients A12 and 
A21 are mechano-chemical coupling coefficients which are 
responsible for energy transduction. 

Looking at the symmetry of the problem, we find that 
V and /cxt transform like vectors for x —x while r and 
A/i are scalars which do not change under inversions. As 
a consequence, the coefficients An and A22 transform as 
scalars while A12 and A21 are vector coefficients. The 
latter can be nonzero only if the system has a polar sym- 
metry. Thus, the polarity of the system (polar filaments) 
is essential for motion to exist. 

Calculating the dissipation rate 11 in the linear regime, 
we find that 11 is positive exactly if the diagonal elements 
are positive, Xu > and if the determinant is positive 



A11A22 — A12A21 > 



(9) 



On general grounds, we expect a symmetry relation be- 
tween the Onsager coefficients if microscopic reversibility 
is obeyed: 



A12 — A21 



(10) 



This is a general result of non-equilibrium thermodynam- 
ics. 



C. Modes of operation 

Different modes of operation of the motor can be dis- 
tinguished by looking at the input and output of energy 
of the system. The dissipation rate 11 corresponds to the 
total fiux of energy to the thermal bath at temperature 
T. Passive regimes of the motor are those cases where 
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both rA/i and fextv are positive: Work performed on the 
system is dissipated and lost. 

More interesting are the active regimes where the mo- 
tor transforms chemical energy into mechanical work or 
vice versa while dissipating only a part of the energy in- 
put. Four such active regimes exist, see Fig. |l]: 

A: rA/i > 0, fc^tv < 0, The motor uses the chemical 
energy of the ATP in excess as input and performs 
mechanical work moving with v > against a neg- 
ative force /cxt < 0. 

B: rA/i < 0, fcxtv > 0, The motor produces ATP, 
although already in excess, from mechanical input 
due to a negative force /oxt < inducing a negative 
velocity v <0. 

C: rA/i > 0, fcxtv < 0, The motor uses ADP in excess 
to perform mechanical work. 

D: rA/i < 0, fcxtv > 0, The motor produces ADP 
already in excess from mechanical work. 



Tjmax 




Figure 1. Operation diagram for an isothermal motor in 
the linear response regime as a function of external force /ext 
and chemical potential difference Afi. General case with four 
different regimes A-D, separated by lines v — and r = 
where the velocity and the fuel consumption vanish, respec- 
tively. The maximal efficiency occurs along a line j^max. 

The different regimes are separated by the lines /ext = 
0, A/x = 0, V = and r = 0. For regimes A and C, where 
the motor performs mechanical work the mechanical ef- 
ficiency is the one defined in Eq. (^): ^7 = — /cxtf /rA/i. 
Similarly, in regimes B and D, where the system performs 
chemical work, the chemical efficiency rj^ = —rAfi/fcxtV 
is more useful. 

Within the linear response regime, the efficiency can 
be calculated using the Onsager coefficients 



where a = /cxt/A/i. If we choose a constant A/i > 0, the 
efficiency vanishes for /oxt — (no work is performed). 77 
becomes positive for /cxt < (note the minus sign which 
indicates that the force is applied in the direction oppos- 
ing movement), reaches a maximum for a certain value 
of the force and becomes zero again at the stall force for 
which V — 0. According to Eq. (pd|), the efficiency is con- 
stant along straight lines /cxt = ciAfi which correspond 
to constant a. Thus, at the origin of the (/cxt, A/i)-plane 
which corresponds to thermal equilibrium and reversible, 
quasistatic operation, the efficiency rj has a singularity 
and is multi-valued. 

Maximal efficiency occurs for a certain value a for 
which drj/da^O. It is given by [|4|j|| 



^max 



(1- VT^)Va 



(12) 



A2ia + A22 



(11) 



Here, A = A^2/(Aii A22). It varies between Tymax — for 
A12 = and ?7max = 1 if A^2 — A11A22. Larger values 
A > 1 violate thermodynamics according to Eq. (^ and 
the Onsager relation (see Eq. (p^). 

These arguments demonstrate that the efficiency van- 
ishes under stalling conditions v = 0. This is an im- 
portant difference from Carnot engine for which the effi- 
ciency is optimized under quasistatic conditions without 
net motion. It results from the fact that the energy trans- 
duction driven by a chemical reaction considered here will 
in general still have a nonzero consumption rate r even 
when motion stops, or in other words v = and r — 
do not occur for the same conditions. 

There is however one limiting case where this is no 
longer true: If A 1, the two lines r — and v ~ 
in the (/ext, A/i)-planc tend towards each other. In this 
limit the chemical reaction and motion are strictly cou- 
pled (i.e. one can't occur without the other) and the 
efficiency reaches the maximum 77 = 1. This situation 
is an idealized case which applies to good approximation 
to polymerization forces and motion generated by poly- 
merization processes as in the case of RNA polymerase 
11,0 



III. TWO STATE MODEL 
A. Transport equations 

We study energy transduction and efficiencies of 
isothermal motors using simple two-state models. The 
motor is characterized by its position x along a one di- 
mensional coordinate describing the polar and periodic 
track. We assume that the motor exists in two different 
conformations or states a = 1,2. The interaction between 
motor and track depends on a and is described by po- 
tentials Wa-{x) with polar symmetry which are periodic 
with period /. 

The role of the chemical reaction is to trigger transi- 
tions between the two states. We introduce the position 
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dependent rate constants uji{x) and (jJ2{x) which char- 
acterize the probabihty per unit time for the transitions 
1 — > 2 and 2 — > 1 at position x, respectively. The proba- 
bihty densities Pi{x,t) and P2{x,t) for the system to be 
at time t at position x in one of the two states obey the 
Fokker-Planck Equations |0] 



dtPl + d^Ji = ~-Ui{x)Pi + LU2{x)P2 
dtP2 + dxJ2 = ^^l{x)Pl - i^2{x)P2 



The particle currents are given by 

J,=r'[~Td,Pa~Pad^W^ 



(13) 



(14) 



where is an effective mobility, the temperature T is 
measured in units of fcs and /cxt is the external force 
introduced above. 

For given rates uJa the system relaxes to a steady state 
with dtPcr = 0. The normalized distributions which sat- 
isfy periodic boundary conditions (Jg dx {Pi + P2) = 1, 
P^(0) = P^{1) and d^PaiO) = d^P^il)) in the steady 
state allow us to calculate the average velocity 



dx (Ji + J2) . (15) 



B. Coupling to a chemical reaction 

We now consider the situation where the transitions 
between states 1 and 2 occur as a result of a chemical 
reaction scheme which we model separately. In order to 
be general and to capture different situations, we consider 
the following scheme: 

ai 

ATP + Ml ^ M2 + ADP + P (16) 
0:2 

71 

ADP + P + Mi ^ M2 + ATP (17) 
72 

Pi 

Ml ^ M2 , (18) 
P2 

where a^, 7^ and (3i denote the forward and backward 
rates, respectively. The reaction pathway a involves 
ATP hydrolysis with chemical free energy gain Afi when 
changing from state 1 to state 2, while pathway 7 in- 
volves hydrolysis in the opposite direction. The transi- 
tions (3 are do not involve a chemical potential difference. 
Chemical kinetics requires 



— = ^{Wi^W2+Ap.)/T 

a2 

71 ^ ^iWl-W2-A^L)/T 

72 

P]_ _ JWi'W2)/T 
132 ~ 



(19) 
(20) 

(21) 



The transition rates can therefore be written as 

=a2e^^^"^"+'^''^/'^+72e(^i"^^"'^'')/^-H/32e(^i"^"^/'^ 

UJ2 =a2+j2+p2 , (22) 

where unknown (/-periodic) functions a2{x), 72(2;) and 
f32 (x) define the conformation dependence of transitions 
rates With these expressions, the net steady state 
ATP consumption rate is given by 



dx [(aiix) - ji{x))Pi{x) - {a2ix) - ^2ix))P2{x)] 

(23) 



C. Detailed balance 

If A/i — 0, the chemical reaction is in equilibrium 
and the transition rates are just thermal fluctuations and 
obey the relation of detailed balance uji/uj2 = exp{{Wi — 
W2)/T). Breaking of detailed balance for A/i 7^ is a 
requirement for spontaneous motion and force generation 
to be possible. In order to quantify the departure from 
thermal equilibrium and the extend to which detailed 
balance is broken, we define the quantity 



il{x) = 1^1 (.t) — 0^2(2;) exp — 



AW{x) 



T 



(24) 



with AW{x) = 1^2(0;) — Wi{x). Detailed balance is 
obeyed only if n{x) = for all x. Using the transition 
rates as given by Eq. (p2[), we find 



n{x) 



'AW/T 



a2 



(x) (e^'^/^-l)+72(2;) (e-^'^/^-1 

(25) 



If Afi 7^ 0, we distinguish two interesting limits: for 
smah A^/T < 1 



ri(x) ~ (02(3;) - 72(a;))e" 



-AW/T 



Afi 



(26) 



indicating that is proportional to Afi. If Afi/T is 
large compared to one, O depends only on the ratio 
k = [ATP]/[ADP][P]: 

n{x)c,{a2{x)e''''°/^k + j2ix)e-^'^°/^k-') , (27) 



where Af/' = fi^j^p — fJ^^ADP ~ Mp- Here, we used the 
relation fii = fi^ + nn[z] where \i\ is the concentration 
of species i, and /i° the so called standard chemical po- 
tential. 
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D. Energy conservation and dissipation 



n = Hi + Ha + n„ 



(35) 



The first law of thermodynamics requires that the en- 
ergy flow through the system is conserved as described by 
Eq. (^. This energy conservation can be derived from 
the transport equations. This leads to expressions for the 
local density of energy dissipation which gives interesting 
insights in how energy transduction is occurring. 

We distinguish two types of dissipation rates: (i) the 
dissipation rates Ha , with a — 1,2 corresponding to slid- 
ing within the potential profiles and (ii) the dissipation 
rates Il^t, with /i = a,/3,7 corresponding to transitions 

In addition to the total dissi- 
^, we introduce local dissipation 



between the two states 
pation rates Ha- and 11 

densities Q^ix) and 0^(x) with 11 ~ dxQ{x) 

For a particle sliding in the potential Wa-{x) with a 
steady state distribution Pa{x) p7|,|2^ 



and 



where 



6ct — —JadxHa 



- dx Ja{x)dxHa{x) 
Jo 



Ha{x) = Wa{x) - /extx + Tln(P,(:E)) 



(28) 



(29) 



(30) 



is an enthalpy whose gradient induces the Fokker-Planck 
current J„: 



Ja = -r'p.dxHa 



(31) 



Therefore, is positive definite as expected for a dis- 
sipation rate. Similarly, the dissipation densities corre- 
sponding to chemical transitions are given by 



Qfi^{PlPl-P2P2){Hl-H2) 



and 



dx Ofi{x) 



(32) 



(33) 



for /i = a, (3,^. Again, IIq, Hp and n.y are positive defi- 
nite as required. For a steady state with periodic bound- 
ary conditions, we can partially integrate Eq. ( |29| ) and 
find together with Eq. (|l|) 



ni+n2 = J^dx {HAJl + H2dxJ2) + fe.tv 

= J^dx{Hi-H2){iViPi-iV2P2)+fc: 



(34) 



Using Eqns. (|l|),(|2|) and (|2|), we find that, the total 
dissipation rate 



satisfies Eq. (0) and energy conservation is obeyed. 

For small A/i and small /oxt, the two state model has 
a linear response regime which obeys the general prop- 
erties required by thermodynamics. In particular it can 
be demonstrated that the model satisfies the symmetry 
relation of Eq. (|l^) as we describe in appendix |^. 
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Figure 2. Three choices of potentials Wi and W2 with pe- 
riod I and transition regions indicated in grey. The position 
a of the maximum of Wi characterizes the potential asymme- 
try, U denotes the potential amplitudes, (a) System A with 
potentials shifted by a distance S and offset Uq- Active transi- 
tions a and thermal transitions /3 are localized within regions 
of size d near the potential minima, (b) System B with sym- 
metric states. The potentials are shifted by a distance of 1/2, 
active transitions a and 7 are chosen such that the system is 
symmetric with respect to an exchange of the two states, (c) 
System C with a flat potential W2 , localized active transitions 
a and non-localized thermal transitions /3. 
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IV. EFFICIENCIES CLOSE TO AND FAR FROM 
EQUILIBRIUM 

A. Specific examples 

We have introduced a general framework which aUows 
us to study a large variety of systems which differ in their 
potential shapes and in the transition rates , 72 and /32 . 
We now discuss three particular examples which we have 
chosen as prototypes to illustrate the physics of energy 
transduction. 

System A is a system with two periodic potentials of 
equal amplitude U which are piecewise linear and which 
are shifted with respect to each other by a displacement (5 
as shown schematically in Fig. ^ (a). Furthermore, they 
differ by a constant value Ui^: W2{x) — W\{x — (5) + C/q- 
The potentials are characterized by the asymmetry pa- 
rameter a which denotes the position of the potential 
maximum of W\. We choose a reaction scheme with 
chemically activated transitions 011.2 between the low en- 
ergy state 1 and the high energy state 2, passive transi- 
tion /3i_2 and 71^2 = 0. The chemical cycle corresponds to 
subsequent transitions a and /3 which we choose localized 
within intervals of size d: 

^ I otherwise, ^^^^ 

localized near the minimum of Wi and (32 (x) = a2{x — S) 
localized near the minimum of W2, see Fig. || (a). Here 
we have for simplicity introduced a single parameter w 
which sets the typical time scale of transition rates. The 
transition rates uja- of system A obey 

ojiix) = a2{x)e('^^-'^-+^f^^/^+a2{x - <5)e(^i~^-)/^ 
uj2{x) = a2{x) +a2{x — S) . (37) 

System A is chosen in such a way that diffusion within 
the potentials is not necessary for motion generation and 
each chemical cycle generates with high probability a for- 
ward step along the x-coordinate. 

System B has different symmetry and different topol- 
ogy of the chemical reaction scheme as compared to sys- 
tem A, see Fig || (b). The two potentials are shifted by 
exactly half a potential period 6 = 1/2: W2 = Wi{x—l/2) 
and Uq = 0. This allows us to introduce a new symme- 
try: the system is invariant under a shift a; — s- a; -I- //2 if 
at the same time the states are exchanged: 1^2. This 
situation is realized by choosing transition rates /3i^2 = 
and 71 (x) = a2{x ~ 1/2) where we localize all transitions 
near the potential minima. We can therefore write for 
system B 

u;i{x) = a2(x)e(^i-'^^+'^^)/'^ + a2ix- 1/2) 

UJ2{x) = 0:2(2;) -I- 02(2; - 

with 02 (x) given by Eq. (|3^). Note, that system B 
involves two active chemical steps per potential period. 



However, because of its additional symmetry it is 1/2- 
periodic. Furthermore, all chemical transitions involve 
ATP hydrolysis, there are no passive transitions. 

System C is shown in Fig. || (c). It is a variant 
of model A with a weakly bound state W2{x) = Uq of 
constant energy. As for system A we choose a reaction 
scheme with 71^2 = and localized active transitions near 
the minima using again definition (|3^) . Since the poten- 
tial W2 is structureless, we assume passive transitions to 
be non- localized with (32{x) = oj. Therefore in system C 

ui{x) = a2(x)e(>^i->^^+^*')/^+c^e('^i-'^^) 
U2{x) = a2{x) +UJ . (39) 

In this case motion generation involves a diffusive step in 
state 2 which we expect to reduce the efficiency of energy 
transduction. 

In order to discuss these models, we identify the rel- 
evant dimensionless parameters: the dimensionless posi- 
tion X = x/l, reduced temperature t = T/U, reduced 
potentials Wa = /U — fcxtl/U and reduced transition 
rates 10^ = uja/^- Eqns. ( p^ and ( p^ can for a steady 
state be written as 

- dsitd^Pi + Pld.iWi) = x{-^lPl + ^2P2) 
-dMP2 + P2dsW2) = xi^lPl - ^2P2) ■ (40) 

The dimensionless parameter 



compares two time-scales: (i) the typical chemical time 
uj~^ and (ii) the typical sliding time in the potentials 
^P/U. For X ^ 1 transitions are fast compared to slid- 
ing while for X ^ 1 sliding is fast. The model is fully 
characterized by the dimensionless parameters Xj T/U, 
Afi/U, a/l, d/l, 6/1 and Uo/U. The results discussed in 
the following section are obtained by numerically solving 
Eq. ( pO| ) with periodic boundary conditions for the three 
different systems. 

B. EflSciencies close to equilibrium 

Numerical examples for the maximal efficiency in the 
linear response regime as a function of temperature are 
displayed in Fig. □ for systems A, B and C and different 
values of X = lo^I /U. They have been obtained by first 
calculating Onsager coefficients from steady state solu- 
tions for small A/i and small /cxt and using Eq. (|T^). 
The orders of magnitude of the efficiency differ for sys- 
tems A, B and C. The efficiency t] depends on x and 
increases in general with increasing x- System B has the 
largest efficiency which approaches 77 ~ 1 for small T/U 
and decreases monotonically as a function of tempera- 
ture. For systems A and C the efficiency has a maximum 
as a function of temperature and vanishes in the limit of 
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small T /U . This indicates in these cases the importance 
of thermal fluctuations for energy transduction. Note, 
that the limit of small temperatures is subtle since in lin- 
ear response A/^ <C T must be obeyed. Therefore, this 
limit corresponds to first sending A/i to zero and T af- 
terwards. Even for small temperatures the system thus 
remains in a regime where thermally activated passage 
over energy barriers rests important. While system A 
can have significant efficiencies of the order of rj ~ 0.06 
in the linear response regime, the efficiency of system C 
which relies on diffusive steps is small (77 ~ 10"''), see 
Fig. I 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 
T/U 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 
T/U 




Figure 3. Maximal efficiency T^max in the linear response 
regime as function of reduced temperature T/U for systems 
A, B and C with a/l — 0.1 as shown in Fig. (a) System 
A with 5/1 = 0.65, Uo/U = 0.4, at different x- (b) Same 
diagram for system B at different x- (c) Same diagram for 
system C with Uo/U = 1.2 at different x- 
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Figure 4. Energy transduction of system (A) with 
a/l = 0.1 and Uo/U — 0.4. (a) Efficiency 77 as a function 
of the external force /cxt for Afi/U = 0.8, x = £.i^^^/U = 200, 
and different reduced temperatures T/U. Broken lines rep- 
resent the approximation discussed in the text, (b) Maximal 
efficiency jymax as a function of A/i/C/ for x = 200 and dif- 
ferent temperatures, (c) Maximal efficiency as a function of 
X for A/^/C/ = 0.8 and different temperatures, (d) Relative 
dissipation rates as a function of A/i for the same system: 
shown are the fraction of energy dissipated by potential slid- 
ing 771 -I- r]2 and the fraction dissipated via active transitions 
77a. For details see text. 
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C. Efficiencies far from equilibrium 



potential M^2- 



We have shown that the two state model transduces 
chemical energy into mechanical work in the linear re- 
sponse regime, however with varying efficiencies. In lin- 
ear response the chemical action represents a small bias 
of the dominant thermal fluctuations. We are now com- 
paring these results with the properties of energy trans- 
duction far from equilibrium. 

System A: Fig. ^ (a) displays the efficiency 77 as a 
function of the applied force for the system A as defined 
in Fig. ^ for A/i/T = 8 and different temperatures. The 
efficiency vanishes for f^xt = as well as for the stall 
force for which the velocity vanishes. For an intermedi- 
ate value of the force, the efficiency reaches a maximum. 
This value r^max is displayed in Fig. ^ (b) as a function 
of Afj.. This diagram reveals the main characteristics of 
energy transduction: for small A/i we find again the non- 
vanishing efficiency of the linear regime. The efficiency 
increases as a function of A//, reaches a maximal value 
and decreases for large A/i to zero. For sufflciently large 
values of A/^t the efficiency increases for decreasing tem- 
peratures and reaches in the example shown a value of 
77 ~ 0.4 for T/U ~ 0.1. The results obtained for different 
temperatures intersect for small A/i which corresponds 
to the observation discussed above that the efficiency in 
the linear regime displays a maximum as a function of 
temperature. Fig. ^(c) shows the behavior of ijmax for 
fixed A^/U as a function of % = S^Pus/U over a range of 6 
decades. The efficiency increases monotonically with in- 
creasing X from zero to a plateau value. As an important 
result we find that the largest values of the efficiency for 
the relevant energy scale T/U ~ 0.1 are of the order of 
77 ~ 0.5 and occur for A/i ^ U comparable to the energy 
difference between the two states at the transition and 
thus far from the linear regime. 

The dissipation rate 11 can, according to Eq. be 
divided into separate contributions of potential sliding 
IIo- and chemical transitions H^, with ^ — a,/3, 7. It is 
useful to define relative dissipation rates 770- = Ilcr/rAfi 
and rja = Ha/rA/i which are analog to the efficiency and 
describe the fraction of dissipated energy relative to the 
consumed chemical work. Note that ij + rji + ri2 + ria + 
rjp -|- 77^ = 1 follows from energy conservation. Fig. ^ (d) 
shows the dominant relative dissipation rates together 
with the efficiency 77. The dominant dissipation is Hi -|- 
112 resulting from potential friction, dissipation IIq of 
chemical transitions, plays a minor role. The dissipation 
77^ corresponding to passive transitions is smaller than 
0.01 and can be neglected. It is therefore not shown. 
The maximum of 77 corresponds to a minimum of 771 -I- 772- 

The main energy loss results from thermally activated 
backward steps. This idea can be directly tested by cal- 
culating the local dissipation density 8i(a;) -I- 02(3;) as 
defined in Eq. (^). This quantity is displayed in Fig. 
^. The plot reveals that maximal dissipation occurs for 
5 — a < X < 5, i.e. along the steep potential slope of the 



1(01+ 02) 




0.2 0.4 0.6 0.8 1.0 
x/l 



Figure 5. Density Qi{x) + 62(2:) as defined in Eq. ( p8[ ) of 
the dissipation rate as a function of position x for system A 
with (////cxt = -0.48, X = 200, A/i/17 = 0.80. 

A second maximum of local dissipation exists for < 
X < a along the steep slope of Wi. In contrast, minimal 
dissipation occurs near the potential minima where tran- 
sitions between states take place. The steep potential 
slopes where the density of energy dissipation is large in- 
deed are accessed via thermally activated backward hop- 
ping events. The probability of such events increases in 
the presence of an "adverse" external force which limits 
the efficiency of the system. 

System B: Fig. ^ displays the same information for 
system B. The diagram reveals that efficiencies are in 
general larger than for system A, reaching values up to 
77 ~ 0.7 for T/U — 0.05. Furthermore, the maximum 
of the efficiency as a function of A/i is less pronounced 
and shifted to small values of A/i as compared to system 

A. One might expect that the dissipation due to passive 
transitions 77^ in system A which does not exist in system 
B could play a role in improving the efficiency of system 

B. However as discussed above rjp can be neglected and 
is thus not responsible for this effect. The main reason 
for the improved efficiency of system B is the fact that 
the effective energy barrier for thermally activated pas- 
sage over the potential maxima is larger in system B as 
compared to system A. Therefore, fluctuations leading 
to "backward steps" in the opposite direction of aver- 
age motion which completely dissipate a consumed ATP 
molecule are less likely. Each active chemical transition 
is thus transduced into work with high probability. Fig. 
^ (c) shows qualitatively the same behavior of the effi- 
ciency as a function of x for system B as compared to 
system A. Also as discussed for system A, the dominant 
dissipation process corresponds to sliding in the poten- 
tials, see Fig.^ (d). Note that the efficiency is larger than 
in system A, which correlates with the increased barrier 
height reducing the probability of backward steps. 

System C: Energy transduction of system C which 
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involves diffusive steps and non-localized de-excitations. 



0.3 



0.1 



/'-^U=0.05 


(a) 


T/U =0.1 0\ 




/ T/U =0.15^^^ 




/ /T/U=a'2&^ 





-0.6 -0.4 
fextl/U 



-0.2 




0.2 0.4 0.6 0.8 
A)I/U 



1.0 1.2 




0.4 0.6 
A|I/U 

Figure 6. Energy transduction of system B with a/l — 0.1. 
(a) Efficiency 77 as a function of the external force /cxt for 
A/i/U = 0.4, X = il'^i^/U = 400, and diflterent reduced tem- 
peratures T/U . (b) Maximal efficiency j^max as a function of 
Ajj,/U for X = 400 and different temperatures, (c) Maximal 
efficiency as a function of x for An/U = 0.4 and different tem- 
peratures, (d) Relative dissipation rates as a function of A^: 
the fraction of energy dissipated by potential sliding rji + 772 
and the fraction dissipated via active transitions rjc, +11^. 

Maximal efficiencies are of the order of 0.02 and thus 



much smaller than those for systems A and B, see Fig. |^. 
As in system A the largest efficiencies occur for A/i ^ T 
and thus far from equilibrium. The reason for the re- 
duced efficiency becomes clear when studying the relative 
dissipation rates shown in Fig. |^ (b): Most energy is in 
this case dissipated by the passive and active transitions, 
potential sliding is less important. In particular, the non- 
localized and passive de-excitations dominate dissipation 
far from equilibrium. 




0.8 1.2 

A|i/u 




0.4 0.8 1.2 
A\x/\j 



T^max 
0.014 

0.01 

0.006 

0.002 
Iff 



T/U =0.10 (C) 



T/U =0.20 



T/U =0.30 
T/U =0.40 




Figure 7. Energy transduction of system C with 
a/l — 0.1. (a) Maximal efficiency as a function of Afi/U 
for X ~ (J?^/U = 5, Uo/U = 1.2 and different reduced tem- 
peratures T/U . (b) Relative dissipation r;i -I- r]2 in the po- 
tentials as well as the dissipation of transitions -qa and 77/3 
corresponding to (a), (c) Maximal efficiency as a function of 
X for A/i/f/ = 1.2 and different temperatures. 

Very striking is the behavior of the efficiency as a func- 
tion of X shown in Fig. ^ (c): The efficiency displays a 
maximum for certain values of % but vanishes both for 
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large and small x- This property reflects the fact that 
a matching of time scales is crucial for this system: The 
life-time in the excited state should be comparable to the 
diffusion-time over a potential period: 

f - T/^uj . (42) 

Therefore, the optimal value of x should behave as Xopt ~ 
T/U which explains the temperature-dependence of the 
maximum in Fig. ^ (c). 

D. Approximation for small forces 

The efficiency far from equilibrium for A/i/T ^ 1 but 
for small forces can be understood by a simple approx- 
imation which we discuss for system A. In the limit of 
large A/i and U /T we ignore spontaneous hopping events 
over the maxima of potential Wi . Every ATP consump- 
tion event corresponds to a transition to the second state 
from which the particle will eventually decay to the first 
state. During this process, it undergoes a forward step 
with probability p+ , a backward step with probability p- 
or it will return to the initial position with probability 
Pq. Here, we have ignored multiple steps, see Fig. ||. 




Here, we have introduced the probabilities Paix) for mo- 
tion in the forward direction after a particle appears in 
state fj at position x. Similarly, 1 — Pa-ix) is the prob- 
ability for backward motion, see Fig. Since the 
two potentials are shifted with respect to each other 
P2{x) — pi{x — S). The probability p+ requires two sub- 
sequent forward movements of this type, p- results from 
two backward movements. As described in appendix 
the probabilities Po-(x) can be calculated approximatively 
for large U/T. Fig. (a) shows the efficiency estimated 
by Eqns. and (Esh together with the numerically ob- 
tained values for comparison. For small forces the agree- 
ment is good, thus confirming our simplified picture of 
energy transduction in this regime. 

V. CONCLUDING REMARKS 

In the previous sections, we have studied the efficiency 
of energy transduction from chemical energy to mechan- 
ical work using a simple two-state model under isother- 
mal conditions. We considered three different examples: 
system A with two shifted potentials and both active 
and thermal transitions between the two states localized 
at the potential minima; system B with an additional 
symmetry between the two states and no passive ther- 
mal transitions; and finally system C with a fiat weakly 
bound state and non- localized passive transitions. We 
demonstrated that energy transduction can be very effi- 
cient in the systems A and B with localized transitions 
and shifted potentials and is at least two orders of mag- 
nitude smaller in system C which requires diffusive steps 
for motion to occur. Interestingly, the largest efficiency 
can occur far from equilibrium. This is in particular the 
case for systems A and C which both are not very efficient 
in the linear response regime. 



Figure 8. Schematic diagram of events after consump- 
tion of one fuel molecule. Forward steps occur with 
probability p+ = pip2, backward steps with probabil- 
ity p_ = (1 — p2)(l — Pi) and neutral steps with 
Po =P2(1 - Pi) +Pi(l -P2) . 

In the presence of an external force /ext i the efficiency 
can thus be estimated as 



fc: 



< a; > 



where 



< X >~ — P-) — v/r 



(43) 



(44) 



is the average displacement per consumed ATP. The 
probabilities p± can be written as 



p+ ^ P2{Q)pi{S) , p_ = (l-p2(0))(l-pi(<5)) 



(45) 



A. Isothermal motors, heat engines and Brownian 
ratchets 

Efficiencies of energy transduction have been studied 
and discussed for a long time. Of particular significance 
is the concept of Carnot which defines the efficiency of 
macroscopic heat engines coupled to two thermal baths 
at temperatures and T+ > T~ as 



^Carnot 



fcxtV 
Q + 



(46) 



where Q+ is the rate of heat transfer from the hot reser- 
voir. This definition then leads to an upper limit of the 
efficiency Tycarnot < {T~^ — T~)/T^ which cannot be sur- 
passed by any heat engine. In order to characterize en- 
ergy transduction in biological systems, a natural choice 
is [ESi 



■n = 



/cxt« 

rA/i 



(47) 
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which we have adopted in this paper, see Eq. (^), and 
which is based on the chemical potential difference be- 
tween fuel and reaction products. As we have discussed, 
this efficiency obeys < 1 in order to satisfy the first 
law of thermodynamics, but there is no nontrivial upper 
bound. In addition to the obvious fact that ?7carnot de- 
scribes a heat engine and rj an isothermal motor, there 
remains a fundamental difference between the two cases: 
the definition of Tycarnot assumes that all heat dissipated 
in the bath T~ is lost. This is true in most practical cases, 
however if the bath at was also used as the hot bath 
of a second heat engine, some of this energy could in prin- 
ciple be reused. Similarly, the definition of 77 takes into 
account the energy of the lower-energy reservoir, thus as- 
suming that the energy of the reaction products remains 
available. One might think that it is possible to avoid 
this difference between the two definitions by choosing: 



(48) 




where would be the chemical potential of the fuel 
(ATP). This definition would share with Carnot's defini- 
tion the viewpoint that the energy of the reaction prod- 
ucts are not useful, and since rj' — rjAii/ iia would lead to 
the upper bound rj' < (pA — ^J'^p) / ^J'A■ Such a choice, how- 
ever, suffers from a serious shortcoming: only chemical 
potential differences are physically meaningful. Depend- 
ing on the state of reference used for defining ^la, the 
value of ^A could be positive, negative or even vanish. 

The example given above demonstrates that compar- 
ing efficiencies can be dangerous as they may be based on 
different definitions corresponding to different points of 
view. This is also the case for ratchet models which have 
been studied in many variants and under widely varying 
physical conditions. All definitions described above have 
been used in the literature: The definition r/camot for sys- 
tems driven by temperature differences 
the definition for rj given by Eq. (^) ||23' 
rj' p^ . Alternative definitions have been proposed for 
situations where the chemical reaction is not fully spec- 
ified ]rT| , p^ . Other definitions of energy transduction 
efficiencies have been used for systems which are driven 
by stochastic or deterministic forces p^ , p^| -pH . Recently, 
Sekimoto has presented a unified picture which includes 
most systems in a common framework However, in 
general, a given definition is adapted to one particular 
physical situation. 



B. The two-state model and biological motors 

One important motivation of this work is to clarify 
the general properties of energy transduction of biolog- 
ical motors. The characteristic behaviors of our system 
A and B with localized transitions and shifted potentials 
are similar to those observed for processive biological mo- 
tors such as e.g. kinesins which move along microtubules 



and for which the consumption of ATP and the subse- 
quent stepping are strongly correlated for small external 
forces [^,^ 32 1 . Kinesin motors consist of two identical 
active head groups which both hydrolize ATP . There 
is evidence suggesting that the motor could "walk" in a 
head-over-head fashion along microtubules, detaching a 
head in the back and reattaching in front of the molecule 
while keeping the second head bound [^,^. In such a 
picture each ATP-hydrolysis cycle leads to a new situa- 
tion where both heads have exchanged their roles and the 
center of mass of the molecule has advanced one filament 
period. This type of motion is captured in a simple way 
in the variant B of our model which is symmetric with 
respect to the two states. Because of this symmetry, both 
states are indistinguishable but the corresponding poten- 
tials are shifted by 1/2: W2{x) = Wi{x -1/2). We there- 
fore identify each of the two states with one kinesin head 
and 1/2 with the filament period, see Fig. ^ (b). Recently, 
system C with a structureless excited state has been used 
for single kinesin heads which were observed to move pro- 
cessively along a microtubule |3^]. Models of this type 
have typically been considered in the context of non- 
processive motors such as myosins which have a weakly 
bound state during their interaction cycle. Myosins inter- 
act with a filament to generate displacements of the order 
of several nm, but they do not continuously move along a 
filament as individual motors since they easily lose their 
track and diffuse away [^Jss). The latter phenomenon 
is not captured in the one-dimensional two-state model, 
however the flat potential of system C requires diffusive 
steps for average motion and the efficiency is therefore 
smaller than for systems A and B. Under physiological 
conditions non-processive motors operate not as isolated 
enzymes but move together in large groups. In this situ- 
ation, however, diffusive steps become unimportant and 
the efficiency becomes large and reaches the same or- 
ders of magnitude as for model A and B described here 

HSjiii. 

When comparing our simple models with biological 
motors, the value of the adimensional parameter x = 
LO^^P /U introduced in Eq. (|4^) is crucial. The relevant 
orders of magnitude for most parameters are well known: 
Energy scales are U ~ lOT typical time scales of 
conformational changes are w ~ 1ms and the relevant 
length scale is / ~ 5 — lOnm ||3^. However, the fric- 
tion coefficient ^ is unknown and difficult to estimate. 
Therefore, we do not know at which value of x biological 
motors operate. The role of x on the functioning of the 
system can be discussed by comparing both the maximal 
efficiency and the dimensionless velocity v/ojl as a func- 
tion of X: see Fig. The diagram reveals that for large 
values of x for which the efficiency is large the veloci- 
ties become small. For small x velocities are optimal but 
efficiency becomes negligible. This observation suggests 
that optimal conditions are obtained in the intermediate 
regime x — 0.1 — 1 where chemical times and sliding times 
along the potential slopes become comparable. If linear 
molecular motors operate in this regime, the microscopic 
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friction coefficient ^ is of the order of 10~^ — 10~^ kg/s. 
If we estimate ^ from simple hydrodynamic arguments 
(^h — 67r77vis07 where ?7vis is some measure of a "local" 
viscosity, we find 77vis — 10— 100 Poise, 10^—10^ times the 
viscosity of water, values compatible with dense macro- 
molecular solutions. Interestingly, this order of magni- 
tude, corresponds to a diffusion coefficient of 4 • 10^ 
m^/s, a value reported recently for single headed kinesin 
p6[ . This observation together with our estimate suggest 
that linear molecular motors are optimized both from the 
velocity and the efficiency standpoint. 
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Figure 9. Maximal efficiencies (broken line) and normal- 
ized velocities (solid line) as a function of x for system B and 
A/i/f/ = 0.6 and T jU = 0.05. 
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APPENDIX A: ONSAGER COEFFICIENTS 

In linear response theory, the behavior of the system 
is completely defined by the Onsager coefficients Ay . 
The Onsager symmetry relation (|lO|) follows from gen- 
eral thermodynamic arguments and the microscopic re- 
versibility. The calculation of Onsager coefficients is dif- 
ficult, however the symmetry relation c an b e verified by 

In section 



general arguments as shown in section A 1 



A 2 we obtain explicit expressions for the coefficients Ay 
for a many-motor system as introduced in IsUll 



1. Symmetry relation for a single motor 

In order to demonstrate the symmetry relation of On- 
sager coefficients for the two-state model, we start from 



the probability distributions at equilibrium (/cxt = 0, 
= 0) as 



with a normalization factor 



„ "'0 



(Al) 



(A2) 



For small fcxtl/T ^ 1 and A/i/T ^ 1 we define devia- 
tions (x) from equilibrium which obey 



P,(x) = 7Ve-^'/^(H-p,(x)) 



(A3) 



Without loss of generality we consider the case where 
only the transitions and Pa- occur but 70- = 0. To 
linear order in A/i, the transition rates can be written as 



(A4) 



with a = a2e-^^/^ and uj{x) = (03 + /32)e-^^/^. Using 
Eq. ( p^ ) we find to linear order 



T 



/ - W'l (x) N 

( e T dxPij + i^{x){pi - P2) = hi{x) 
dx [e ^ dxP2j -uj{x){pi -P2) = h2{x) . (A5) 



The fields hg- (x) are nonzero in the presence of mechani- 
cal or chemical forces: 



hi{x) 



W-i (x) 

-a{x)^^l/T - dxe ^/ext/^ 



h2{x) — a{x)A^/T — dxG 



fcxt/^ 



(A6) 



Eq. (A5) represents a linear relation between p^. and ha 
which can be inverted and which thus defines a response 
kernel 

P<y{^)^y^j dx'xap{x,x')hp{x') . (A7) 
„ Jo 



This allows us to express the velocity and the fuel con- 
sumption rate within linear response theory: 



dx 



dx a{x) [(pi 



W2 



- TdxP2) 

-P2) + Am/T] 



(A8) 
(A9) 



The Onsager coefficients Ay can be written in terms of 
the response functions Xap{x,x'). In particular, we find 
for the coefficients of mechano-chemical coupling 
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Ai2 



dAfi 

dx / dx' 
Jo 



dtPl+vd^Pl = -LUl(x)Pl+LU2{x)P2 
dtP2 + vd,P2 - LJi{x)Pi - L02{x)P2 



(A17) 



+ ^-^e 5^(9:j:X2i(a;,x') - 52;X22(a;,a;')) a(a;') 

(AlO) 



A21 



dr 

dfcKt 

I pi 

dx / dx' 
Jo 



the velocity v is determined by the force balance condi- 
tion 



[ dx (PAWi + P2d^W2) 

Jo 



(A18) 



Steady state distributions Pi and P2 = 1/Z — Pi are so- 
lutions to 



{Xiiix,x') ~ X2i{x,x'))^ ^da:'e T 



+ (Xi2(a;,a;') - X22(a;,x'))^ ^d^-e t a{x) 



(All) 



Performing a partial integration in Eg. (AlO), we find 
that the Onsager symmetry relation ( |lO| ) is satisfied ex- 
actly if the response functions obey the symmetry rela- 
tion 



Xap{x,x') = Xpa{x\x) 



(A12) 



This symmetry relation follows from the hcrmiticity of 
the linear operator C defined in Eq. (A5) which can be 
expressed as 



C 



where 



with 



—LU C2 + 



(A13) 



(A14) 



(A15) 



The operator C is hermitian since the matrix (A_14) is 
symmetric and Ci itself is hermitian. The latter is easily 
verified by partial integration: 



dxq{x){dj:e t dxp{x)) 



dxp{x){dxe T d^q{x)). 



(A16) 



2. Onsager coefficients for many rigidly coupled 
motors 

Onsager coefficients can be calculated explicitly for a 
model of rigidly coupled motors. In this model, the two 
states are defined in the same way as before, but the 
current term in Eq. ( p^ ) is replaced by a convective 
term since all particles move with the same velocity. The 
transport equations are given by |ll|,p[ 



vdxPi — (wi + i^2)Pi +1^2/1 



(A19) 



Using a power expansion in the velocity. Pi can be writ- 
ten to lowest order 



P^{x) = Pf\x) + P[^\x)v + 0{v^) 
with P^^ = LO2/ {loi -\- 102)1, and 



^1^ 



(1) 



1 



bJl + UJ2 



d.Pl 



(0) 



(A20) 



(A21) 



As in the last section we use a reaction scheme with 7^ = 
in order to keep the expressions simple. Also, without 
loss of generality we consider the case where ai depends 
on A/i but a2 remains constant. For small Afi/T <^ 1, 
we express the transition rates of Eq. (p^) as 



u;i=c^2e-^^/^(l 



(A22) 



where a = a2A/z/a'2P- The force-velocity relation for 
small V is given by 



/cxt = /(°) + (r^ + /W)t^ + o(^'^) 

with coefficients 



(A23) 



/(") = - / dxP^^'^dAW . (A24) 
Jo 

which depend on A^. We can now calculate the Onsager 
coefficients. The effective friction 



All 



dv 



Ap,=0 



9/oxt 

can be determined from Eq. ( A23D 
1 



A 
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Similarly, 



leads to 



dv 



IT uj2{l + e-^w/TY 

dfcxt 



dAfj. 



-All 



/oxt=0 



dA^i 



(A25) 



(A26) 



(A27) 



•0=0 
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Ai2 — — Ai] 



a/(o) 



dAfi 



Ati=0 



An f' a2e-^'^/^d,AW 

IT Jo ''''6.2(1 + 6-^^/^)2 

(A28) 



The second cross-coefficient 
dr 



A 



21 



A/j=0 



Or 

All 
ov 



(A29) 



A/j=0 



is determined from the fuel c onsum ption rate r(w,A/i). 
Using Eqns. (psl), (|A29D and (|A2lD, we obtain 



9r 

9u 



= / (q;2 + ai)9i,Pi] 

=0 "'0 



and thus as required A12 = A21. FinaUy, 



A2 



dr 



dAfi 

^12 ^ 

All 9A/i 



A 
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/oxt=0 

dr 



dr 
dv 



dr 



Am=0 



(A31) 



u=0 



which requires to calculate 



dr 



dA^i 



v=0 -JO 



^(ai+a2) + Pl^^(ai+a2) 



(A32) 



Using Eqns ( ^21 ) and (|23|), we obtain 



i9r 



dx 



-\W/T 



(1 + e~^^'T) 



UJ2 



(A33) 



Note that both An and A22 are positive while A12 can 
have either sign. 



APPENDIX B: DIFFUSION CLOSE TO A 
POTENTIAL MAXIMUM 



In Section IVD, we introduced the probabilities Po- (a^) 



that a particle which initially starts at position x close to 
a potential maximum will finally escape in the positive 
direction. We will calculate this probability for a piece- 
wise linear potential as shown in Fig. 10 in the limit 
where the potential slopes extend to infinity which cor- 
responds to large potential amphtudes U/T » 1. We 
consider the Fokker-Planck Equation 



dtP + d^J = 



(Bl) 



with 



J^~r'[Td,P + Pd.W-Pf,^t] , (B2) 

for initial condition 

P{x,t^0)^5{x-xo) . (B3) 

In order to determine what fraction of particles move to 
the right after a long time, we define the Laplace trans- 
form 



P{x,s) = / dtP{x,t)e~ 



(B4) 



of the distribution and J(x, s) = -^-^[Td^P + Pd^W - 
Pfext] of the current. 



P(x,t=0)=5(x-xo) 




Figure 10. Schematic representation of diffusion near a po- 
tential maximum which can be divided in three different re- 
gions A, B and C. For a particle which initially appears at 
a; = xo, we are interested in the probability p that it will fi- 
nally move forward. This probability is related to the current 
Jc which can be calculated in the limit of long times. 

The average number of particles which pass at position 
X after long times is given by 



dtJ{x,t) = J{x,s = 0) 



Noting that 



d 



dte^'^ — Pix^t) = -P(x,0) + sP(x,s) 

we obtain an equation for P: 

\Td,P + {d:,W - foxt)P 



(B5) 



(B6) 



sP 



-Pix,0) 



(B7) 



Since we are interested in s = we have to solve 

9, hd^P + {d^W - Ut)P] = -^Six - xo) . (B8) 
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All quantities of interest can be easily calculated if we we 
assume a piecewise linear potential 



W{x) = 



U + f-X 
U-f+x 



X <0 
x> 



(B9) [2: 



with the potential slopes /_ — U/a and /+ = U/{1 — a). 
We distinguish three different regions A, B and C along 
the a:-axis, see Fig. 10. Within each region, the solution 
to Eq. @ is 

P{x) = Co + Cie-(^(^)-^^-')/^ , (BIO) 

with two constants Cq and Ci which have to be deter- 
mined for each of the three regions. We denote the corre- 
sponding solutions Pa, Pb and Pq. Since we are looking 
for solutions which do not diverge for a; — + ±00, we have 
Ci = Ci = in regions A and C and therefore Pa — Cq 
and Pc — Cq . This boundary condition for large x can 
also be derived more carefully by first imposing the con- 
dition 



lim P{x, s) 

; — >ioc 







(Bll) 



for s > and taking the limit s — s- afterwards. Addi- 
tional boundary conditions are the conditions of continu- 
ity of P{x, 0) at a; = and x = xq 



(B12) 



PAiO) = PB{0) 

pBixa) ^ Pc{xo) , 
and the matching conditions 

d^pAiO) = d^pBiO) - Pa(0)(/+ + f-)/T 
d^pBixo)^d^Pc{xo)+^/T , (B13) 



which follow from integrating Eq. (B8) at the singulari- 
ties of the potential slope at x = and the delta- function 
a.t X = Xq. With these conditions, all free parameters can 
be determined and we obtain 



ck 



./+ + /cxt 



= -(/++/oxt)a:o/T 



/--/■ 
/- + /+ 



The corresponding currents J{x, s - 
region A and C: 

Ja = -£.^^i.f- - fc^t)CQ 

Jc ^ r\f+ + fc.t)c^ 



(B14) 

(B15) 

0) are constant in 

(B16) 
(B17) 



which satisfy the normalization condition Jc ~ Ja — 1- 
The probability p{xo) for forward motion of a particle 
which initially was at xq is equal to Jc'- 



/ \ -, /— /cxt - 

p[xa) = 1 - - — — 
J- + J+ 



The probabilities introduced in Eq. 

P2{0)^p{l-5),p^{5)=p{5-a). 



(B18) 
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